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(N ! Abstract 

According to the driving potential of a fissile system, we propose a phenomenological fission 

potential for a description of the pre-neutron emission mass distributions of neutron-induced ac- 

tinide fission. Based on the nucleus-nucleus potential with the Skyrme energy-density functional, 

the driving potential of the fissile system is studied considering the deformations of nuclei. The 

energy dependence of the potential parameters is investigated based on the experimental data 

for the heights of the peak and valley of the mass distributions. The pre-neutron emission mass 

distributions for reactions 238 U(n, f), 237 Np(n, f), 235 U(n, f), 232 Th(n, f) and 239 Pu(n, f) can be 

reasonably well reproduced. Some predictions for these reactions at unmeasured incident energies 

' are also presented. 
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I. INTRODUCTION 



Nuclear fission is a field of very intense studies in the recent decade lKZj]. One of the 
most interesting characteristics of neutron-induced fission is the huge difference of the mass 
distribution of the fission fragments for different nuclei and the dramatic change of the mass 
distribution with the variation of the incident energies of neutron. It is still far from clear 
how the parent nucleus transforms into a variety of daughter pairs. The highly excited 
primary fission fragments, whose mass distributions are called pre-neutron emission mass 
distributions, are de-excited by the emission of prompt neutrons, and followed by prompt 
7-rays to form the primary fission products. The primary fission products are usually highly 
neutron-rich and unstable, and gradually evolve to the secondary fission products through 
the emission of delayed neutrons and the radioactive /3-decay. We will focus on the study 
of the pre-neutron emission mass distribution of the primary fission fragments in this work. 
The precise calculation of the pre-neutron emission mass distributions is of great importance 
for understanding the fission process and for describing the yields of the fission products. 

The measured mass distribution of the fission fragments can be reasonably well repro- 
duced with some empirical approaches or some systematical methods. Liu et al. had de- 
veloped systematics of mass distributions for neutron-induced 238 U fission [8] and of inde- 
pendent yields for neutron-induced 235 U fission [9|. They also presented the evaluation data 

JjJ and Wahl Q had 



1 Of ] and the adjusted data for several actinides. Katakura [12I, 



fitted the experimental data for actinide nuclei with 3 to 7 Gaussian functions. Kibkalo's 
phenomenological model was designed to study the dependence of the mass distribution on 
;he transferred angular momentum, and was later adapted for predictions of fission yields 
JjJ. A new systematics for fragment mass yields of target nuclei from Th to Bk at incident 
particle energies between 5 and 200 MeV was developed by Gorodisskiy et al. through inde- 
pendent fission modes [3] . The mass distributions predicted with the systematical methods 
mentioned above are generally described by a series of Gaussians, and the model parame- 
ters are obtained through fitting the experimental data. However, available experimental 
data for energy-dependent neutron-induced fission yields, especially the pre-neutron emis- 
sion mass distributions, are not sufficient enough for the development of global systematics, 
which results in great difficulties for the predictions of the mass distributions at unmeasured 
energies. 
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For a microscopic description of the mass distribution of nuclear fission, the precise cal- 
culation of the potential energy surface seems to be required. Unfortunately, the micro- 
scopic calculation of the potential energy surface of a fissile system is very complicated and 
time-consuming. Some phenomenological approaches are still required for the quantitative 
description of the energy dependence of the mass distribution at present. It is known that 
the shell and pairing effects play a key role for the fission and quasi-fission process. It is 
found that the quasi-fission mass distribution in fusion reactions leading to the synthesis 
of super-heavy nuclei can be reasonably well described by the driving potential in the di- 
nuclear system (DNS) model [l?], Q], since the shell effects of reaction system are effectively 
involved in the driving potential via the Q-value of the system. It is therefore interesting to 
investigate the mass distribution of a fissile system based on its driving potential. Because 
the deformation effect of a nuclear system influences the process of fusion and fission jlo| . 
it is expected that the deformations of nuclei play a role for a reliable calculation of the 
corresponding driving potential. 

In this work, we attempt to propose a simplified fission potential with a few well- 
determined parameters for quantitatively describing the pre-neutron emission mass distribu- 
tions of neutron-induced actinide fission, by combining the corresponding driving potential 
of the fissile system. This paper is organized as follows: In Sec. II, we briefly introduce the 
calculation of the driving potential of a fissile system. In Sec. Ill, we introduce the fission 
potential and its parameters. In Sec. IV, the comparisons between the predicted results and 
the measured data of the pre-neutron emission mass distributions for the reactions 238 U(n, 
f), 237 Np(n, f), 235 U(n, f), 232 Th(n, f) and 239 Pu(n, f) are presented. Finally, the summary 
and discussion is given in Sec.V. 

II. DRIVING POTENTIAL OF A FISSILE SYSTEM 

Assuming that a compound nucleus separates into a pair of nuclei in the fission process, 




(1) 



the corresponding Q-value of the system can be expressed as, 




(2) 
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Here, E(Ai, Zi) denotes the energy of a nuclear system with mass A\ and charge Zj. For 
a description of the potential energy surface of a fissile system around scission point, the 
dinuclear system (DNS) concept may be used, omitting the excitation energies of the fission 
fragments for simplicity According to the DNS concept, each fission fragment at the scission 
point retains its individuality in the evolution of the DNS. This is a consequence of the 
influence of the shell structure of the partner fragments since the collective kinetic energy 
of the fission fragments is low around the scission point. Based on the DNS concept, the 
driving potential of a fissile system is expressed as 

D = Q + B (3) 

in this work. Bo denotes the Coulomb barrier height in the interaction potential between the 
fragment pair. In this work, the interaction potential V(R) at a center-to-center distance 
R is calculated by using the Skyrme energy-density functional together with the extended 



Thomas- Fermi (ETF) approximation 201 ] . 



Fig. 1 shows the driving potential for U(n, f). The solid circles denote the sum of the 



Q-value and the average total kinetic energy TKE of fission fragments at an incident energy 
of 2 MeV. One sees that there exists a valley at ~ 140 for the mass number of the heavy 
fission fragments. The open circles denote calculated driving potential with the deforma- 
tions of nuclei being taken into account. Here, the deformations of nuclei at their ground 
state are taken from the calculations of the finite range droplet model [21 1. We consider 
the tip-tip orientation and simultaneously consider the dynamical octupole deformation of 
fragments (empirically set \^\ = 0.08) in the fission process around the scission point, since 
one obtains a lower Coulomb barrier at the tip-tip orientation. We note that the measured 



data for Q + TKE are roughly reproduced by the calculated driving potential. The crosses 
denote the results without the deformations of nuclei being taken into account. From the 
comparison, one learns that the deformations of fragments play an important role for a rea- 
sonable description of the fission mass distribution and the total kinetic energy of fragments. 
The solid curve denotes the result from an empirical fission potential which will be discussed 
in next section. 
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FIG. 1: (Color Online) Driving potential for 235 U(n, f). The solid circles denote the sum of the 
Q-value and the average total kinetic energy TKE of fission fragments. The crosses and open 
circles denote the driving potential D without and with the deformations of nuclei being taken 
into account, respectively. The solid curve denotes the results of kU + 10 with k = 1 MeV. U{A) 
denotes the empirical fission potential which will be discussed in next section. The sub-figure shows 
the corresponding density contour plots of a typical primary fragment pair 96 Sr+ 140 Xe. 

III. FISSION POTENTIAL AND ITS PARAMETERS 

For a more quantitative description of the mass distribution of primary fission fragments, 
we further propose an empirical fission potential considering the calculated driving potential. 
One expects that the pre-neutron emission mass distribution, i.e., the mass dependence of 
the primary fission fragments are strongly dependent on a corresponding fission potential. 
We assume that the pre-neutron emission mass distributions of low-energy neutron-induced 
actinide fission can be approximated described by using a simplified fission potential U(A), 

P{A) = Cexp[-U{A)]. (4) 

Where C is the normalization constant, and the variable A denotes the mass number of 
the primary fragment. Considering the double-humped mass distributions of low-energy 
neutron-induced actinide fission, we describe the phenomenological fission potential U(A) 
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by using three harmonic-oscillator functions, i.e., 

u^A-Ai) 2 A<a 
U(A) = < —uq{A — A ) 2 + R a < A < b (5) 
u 2 {A-A 2 ) 2 A>b. 

Where, Ai, and A 2 are the positions for the peaks of the light and heavy fragments of the 
pre-neutron emission mass distributions, respectively. According to the calculated driving 
potentials, we find that the valley for the mass number of heavy fragment locates A 2 ~ 140 
for neutron-induced actinide fission in general. Therefore, we set A 2 = 140 in the calculation 
for simplicity. We have checked that the calculated mass distribution of fission fragments 
do not change appreciately if the value of A 2 is slightly changed. A = A CN /2 denotes the 
corresponding position for symmetric fission. Here, Acn is the mass number of the fissile 
nucleus. Considering that the fission potential is a smooth function, the coefficients in Eq. (5) 
can be derived as 
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Uq = 
Ul = 

U2 = (A 2 -A )(A 2 -bY (6) 

with A\ = Aqn — A 2 and b = - — - — — — — — + A . The potential parameters a and R 

A 2 — Aq 

will be discussed later. 

The total mass distributions of the binary fission fragments should be normalized to 
200%. The normalization constant C can therefore be analytically expressed as 

200% 200% 
1 = / °° exp[-U(A)]dA ~ Io + h + lY U 



with 



Io = ^^{erfi[(Ao-a)V^]+erfi[(6-A )V^]}, 



h = -^{erf^a-AOV^il+erf^V^i]}, 



h = ^{l + erf[(A 2 -6)V^]}- (8) 
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Where erf(x) and erfi(x) denotes the error function and imaginary error function, respec- 
tively. We also assume that P(Ai) = P{A-i), i.e., for the pre-neutron emission mass distri- 
butions the height of the peak of the light fragments equals to that of the heavy fragments. 
The parameter a can be uniquely determined by the normalization constant C and P(A 1 ). 




FIG. 2: (Color Online) Values of peak P{A\) and valley P(Aq) of the pre-neutron emission mass 
distributions for reaction 238 U(n, f) as a function of incident energy of neutron. The data are taken 



from Ref. 



221. The solid lines denote the results in this work. 



The potential parameter R is defined as 

*-*m 

We find that the heights of the peak and valley of the mass distributions change linearly 
with the low incident energies of neutron in general. The energy dependence of P(Aq) and 
P(Ai) are written as, 

P{Aq) = a + (3 E n , 

P(A X ) = a x + PtEn. (10) 

Here, E n denotes the incident energy of neutron. The parameters atQ, (3 , a,\ and /3i are finally 
determined by the experimental data for P(Aq) and P(A\). The potential parameters for 
different reaction systems are listed in Table I. 



TABLE I: The potential parameters adopted in this work. 



reaction 


a i 


ft 


Q'O 


A> 


232 Th(n, f) 


7.0900 


-0.0631 


-0.0656 


0.0350 


235 U(n, f) 


6.6550 


-0.0952 


-0.0042 


0.0237 


237 Np(n, f) 


6.4132 


-0.1008 


-0.0113 


0.0338 


238 U(n, f) 


6.5508 


-0.1090 


-0.0208 


0.0196 


239 Pu(n, f) 


6.1293 


-0.1517 


0.0067 


0.0343 




FIG. 3: (Color Online) Pre-neutron emission mass distributions at an incident energy E n = 0.5 
MeV for reaction 235 U(n, f). The scattered symbols denote the experimental data which are taken 



from Ref. 



23J. The solid curve and bars denote the calculated results with the empirical fission 



potential U (A) and the driving potential, respectively. 



IV. RESULTS 



In this work, we first investigate the energy dependence of the potential parameters. In 
Fig. 2, we show the values of P{A{) and P(Aq) in the pre-neutron emission mass distribu- 
tions for the reaction 238 U(n, f) as a function of the incident energy of neutron. One can 
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FIG. 4: (Color Online) Pre-neutron emission mass distributions at incident energies from 1.3 MeV 
to 6 MeV for reaction 238 U(n, f). The scattered symbols denote the experimental data which are 
taken from Ref. [zj (squares) and Ref. [2J] (circles), respectively. The solid curves denote the 
calculated results in this work, (d), (h) and (1) show the predicted results at three unmeasured 
energies 2.2, 4.2 and 6.0 MeV, respectively. 

see that the values of P(Ai) and P{A ) linearly change with the incident energy in general, 
which could provide us with some useful information for calculating the pre-neutron emis- 
sion mass distributions at unmeasured energies of E n = 1 ~ 6 MeV. From the potential 
parameters listed in Table I, one can see that the values of ol\ and f3\ decrease with the mass 
number of the fissile nuclei in general. We also note that the potential parameters, such as 
A 2 , P(Ai) and P(A ), are different in different models and the energy dependence of A 2 is 
weak. 

In Fig. 3, we show the comparison of the calculated pre-neutron emission mass distribu- 
tions P{A) for 235 U(n, f) with the the empirical fission potential and the driving potential, 
respectively. The solid curve denotes the result with U(A) in Eq.(4). The bars denote the 
corresponding result with the obtained driving potential in Fig. 1. Here, the mass distri- 
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FIG. 5: (Color Online) The same as Fig. [H but for reaction Np(n, f) at incident energies from 



0.3 to 6.0 MeV. The experimental data (square dot) are taken from Ref. 



25] 



bution is roughly estimated by using a formula P(A) oc exp(— D/k) based on the obtained 
driving potential D with k — 1 MeV and considering the normalization. One sees the the 
positions and widths of the peaks for mass distribution can be reasonably well reproduced 
with the driving potential. The large fluctuation is due to that the temperature dependence 
of nuclear structure effect is not considered yet. With the empirical fission potential, the 
description of the measured mass distribution for the primary fission fragments in 235 U(n, 
f) can be significantly improved. 

In Fig. 4, we show the pre-neutron emission mass distributions P(A) at incident energies 
from about 1 to 6 MeV for the reaction 238 U(n, f). The scattered symbols denote the 
experimental data which are taken from Ref. 22j (squares) and from Ref. 24j (circles), 
respectively. The solid curves denote the calculated results in this work. The potential 
parameters «o, Po, «i and j5\ adopted in the calculations are listed in Table I. Fig. 5, Fig. 6 
and Fig. 7 show the pre-neutron emission mass distributions P(A) for the reactions 237 Np(n, 
f), 232 Th(n, f) and 239 Pu(n, f), respectively. One can see that the experimental data can 
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FIG. 6: (Color Online) The same as Fig. HJ but for reaction U(n, f) at incident energies from 



0.5 to 5.5 MeV. The experimental data (square dot) are taken from Ref. 



23] 



be reproduced reasonably well, which indicates that the fission potential proposed in this 
work is reasonable. For some unmeasured energies, we also present the predictions from this 
approach. 

V. SUMMARY AND DISCUSSION 

In this work, we proposed a phenomenological fission potential based on the corresponding 
driving potential for quantitatively describing the pre-neutron emission mass distributions of 
neutron-induced actinide fission at incident energies of neutron at a few MeV. Based on the 
nucleus-nucleus potential with the Skyrme energy-density functional, the driving potential 
of the fissile system is studied considering the deformations of nuclei. The measured data for 



the sum of the Q-value and the average total kinetic energy TKE of fission fragments at an 
incident energy of 2 MeV in 235 U(n, f) can be reasonably well reproduced by the calculated 



II 




FIG. 7: (Color Online) The same as Fig. [H but (the upper panel) for reaction 232 Th(n, f) and 
(the lower panel) for reaction 239 Pu(n, f) at different incident energies. The experimental data are 
taken from Refs. 



driving potential. We also learn that the deformations of nuclei play an important role for a 
reliable calculation of the driving potential. With a systematic study on the reactions 238 U(n, 
f), 237 Np(n, f), 235 U(n, f), 232 Th(n, f) and 239 Pu(n, f), we find that the experimental data 
of these reactions can be reproduced reasonably well with the proposed fission potential. 
This investigation is helpful for further describing the yields of the fission products. By 
combining the radial basis function approach [28| , the accuracy and predictive power of the 
model could be significantly improved. In addition, a more microscopic description of the 
potential parameters and the temperature dependence of the driving potential should be 
further investigated. The study on these aspects is under way. 
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